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We study in detail the Quantum Field Theory of mixing among three generations of Dirac 
fermions (neutrinos). We construct the Hilbert space for the flavor fields and determine the gener- 
ators of the mixing transformations. By use of these generators, we recover all the known param- 
eterizations of the three-flavor mixing matrix and we find a number of new ones. The algebra of 
the currents associated with the mixing transformations is shown to be a deformed su{3) algebra, 
when CP violating phases are present. We then derive the exact oscillation formulas, exhibiting new 
', features with respect to the usual ones. CP and T violation are also discussed. 
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I. INTRODUCTION 

In rece nt y ears, there has been much progress in the understanding of flavor mixing in Quantum Field Theory 
(OFT) . The original discovery of the unitary inequivalence of the mass and the flavor representations in QFT 

J> 0, has prompted further investigations on fermion mixing as well as on boson mixing It has emerged 

that the rich non-pcrturbative vacuum structure associated with field mixing leads to relevant modification of the 
00 fiavor oscillation formulas, exhibiting new features with respect to the usual quantum-mechanical ones ]T8[ |. Some 
topologically nontrivial features, such as the occurrence of a geometric (Berry- Anandan) phase in field mixing has 
">! been also pointed out | p^ . 

In this paper we study in detail the case of three flavor fermion (neutrino) mixing. This is not a simple extension 
(—V of the previous results since the existence of a CP violating phase in the parameterization of the three-flavor 

mixing matrix introduces novel features which are absent in the two-flavor case. We determine the generators of 
' the mixing transformations and by use of them, we recover the known parameterizations of the three-flavor mixing 
' , matrix and flnd a number of new ones. We construct the flavor Hilbert space, for which the ground state (flavor 
vacuum) turns out to be a generalized coherent state. We also study the algebraic structure of currents and charges 
(-H ' associated with the mixing transformations and we find that, in presence of CP violation, it is that of a deformed 
• • ■ su(3). The construction of the flavor Hilbert space is an essential step in the derivation of exact oscillation formulas, 
which account for CP violation and reduce to the corresponding quantum-mechanical ones in the relativistic limit. 
The paper is organized as follows. In §11 we construct the Hilbert space for three-flavor mixed fermions. In §111 we 
$H ' study the various parameterization of the unitary 3x3 mixing matrix obtained by use of the algebraic generators. In 
5^ ! §IV, we study the currents and charges for three-flavor mixing, which are then used in §V to derive the exact neutrino 
oscillation formulas. Finally in §VI, CP an T violation in QFT neutrino oscillations are discussed. §VII is devoted to 
conclusions. In the Appendices we put some useful formulas and a discussion of the arbitrary mass parameterization 
in the expansion of flavor flelds as recently reported in |j7|,p| . 
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II. THREE FLAVOR FERMION MIXING 

We start by considering the following Lagrangian density describing three Dirac fields with a mixed mass term: 

C{x) ^^>fix){i^-M)^iJfix), (1) 

where = (z/g, j/^, 1/^) and M = MMs the mixed mass matrix. 

Among the various possible parameterizations of the mixing matrix for three fields, we choose to work with the 
following one since it is the familiar paramctrization of the CKM matrix ||2^,|2^ : 

C12C13 S12C13 siae""*^ 

*/(a;) =Z^*„i(a;)= ( -S12C23 - ci2S23Si3e** C12C23 - si2S23Si3e*'' S23C13 | *m(a;), (2) 

S12S23 - Ci2C23Sl3e*'' -C12S23 - Sl2C23Sl3e*'^ C23C13 
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with Cij = cos 9ij and = sin^y , being 9ij the mixing angle between z/^, Vj and Vf^^ = (i^i, 1^2, t's)- 

Using Eq. (^) , we diagonalize the quadratic form of Eq. (p , which then reduces to the Lagrangian for three Dirac 
fields, with masses mi, m2 and 7713: 



where = dmgfmi, TO2, ms). 

Following Ref . 0] , we construct the generator for the mixing transformation (||) and define]^ 



where (cr, i) = (e. 1), (/i, 2), (r, 3), and 



where 



Gi2(0 =exp [^?12il2W 

G23(i) = exp 6*23^23(0 
Gi3W=exp[0i3Li3(<5,t) 



G(,(t) =G23WGi3(0G'i2W, 

; Li2{t) = I (fx v\{x)v2{x) - h'\{x)vi{x) 



L2z{t) 



d-^x 



v\{x)vz{x) - vl{x)v2{x) 



Li-i{S,t)^ / d^x i^l{x)i^3{x)e - vl{x)vi{x)e 



iS 



(3) 

(4) 
(5) 

(6) 
(7) 
(8) 



It is evident from the above form of the generators, that the phase S is unavoidable for three field mixing, while it 
can be incorporated in the definition of the fields in the two flavor case. 

The free fields I'i (i==l,2,3) can be quantized in the usual way [p2| (we use t = xq): 



^«(^) =Y. +vL^M0 



1,2,3, 



(9) 



with ul^{t) = e-*'"'=.'*Mi, ^, vl^^{t) = e*"'=.'*wi;,^ and w^,, = ^f^T' 



The vacuum for the mass eigenstates is 



denoted by \0)m'- Q^kJO)™ = /5ki|0)m = 0. The anticommutation relations are the usual ones; the wave function 
orthonormality and completeness relations are those of Ref. ^ . 

The main result of Ref. is the unitary inequivalence (in the infinite volume limit) of the vacua for the fiavor 
fields and for the fields with definite masses. There such an inequivalence was proved for the case of two generations; 
subsequently, in Ref. | ^ , a rigorous general proof of such inequivalence for any number of generations has been given 
(see also Ref. [|3[). Thus we do not need to repeat here such a proof and we define the the flavor vacuum as: 



|0(<))/ ^ G,\t) |0). 



(10) 



The form of this state is considerably more complicated of the one for two generations. When (5 = 0, the generator 
Go is an element of the S'f/(3) group (see §IV) and the flavor vacuum is classified as an SU (3) generalized coherent state 
a la Perelomov p3[ |. A nonzero CP violating phase introduces an interesting modification of the algebra associated 
with the mixing transformations Eq.(H): we discuss this in §IV. 

By use of Gg{t), the flavor fields can be expanded as: 



'^Let us consider for example the generation of the first row of the mixing matrix U. We have (see also Appendix B) 
dve/de23 = 0; and 9!/e/aei3 = G-^^G-^^ui, Li3]Gi3Gi2 = G-^^G^^^e-'^usGisGii, thus: 

d^Ve/dOlz = -Ve => fe = f (Ou) COS 0vi + g{9x2) 0-iz\ 

with the initial conditions (from Eq.(^): /(S12) ~ i^e|si3=o and (7(^12) = du^/ d9ri\ei3=o ~ e~'*i^3 . We also have 

f {612)/ dels = -f{0i2) => f {612) = A cos 012 +B sin 012 

with the initial conditions A = Ve\e=o ~ i^i and B = df(0i2)/d0i2\e=o ~ 1^2, and — (S12, ^13, ^23). 
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^-(^) = E / [K/tXAt) + V-K^it)f^-iJ^)] ^"^ '^ ' - (e, 1), (m, 2), (r, 3) . (11) 

r •' 

The flavor annihilation operators are defined as a'^^{t) = G^"^ (i)a{^ jG'e(t) and P^]^^{t) = Gg^{t)f3'^^^ fieit)- They 
clearly act as annihilators for the fiavor vacuum Eq.([l0|). For further reference, it is useful to list explicitly the flavor 
annihilation/creation operators (see also Ref. [|l|). In the reference frame k = (0,0, |k|) the spins decouple and their 
form is particularly simple: 

«i;.e(0 = C12C13 + ([/^^(O a^^^ + e'-Vt^it) p^,) + e~^' s,, {u^;{t) a^^^ + ,^Vt^{t) (i^^) , (12) 

aLp(^) = (ci2C23 - e'^ S12S23S13) ai,2 - (si2C23 + e'^ C12S23S13) (u^2{t) "Li " ^''"^12(0 P'iks) 

+ S23C13 (U^s*it) <3 + ^'■^23(0 ' (13) 

<,r{t) = C23C13 aL.S - (C12S23 + S12C23S13) (u^sit) aj, 2 - e'^^2''3(0 /3-k,2) 

+ (S12S23 - e^' C12C23S13) (U^t) al, - e^V^.it) f^^i) ' (14) 



P-uJt) = (C12C23 - e-^' S12S23S13) /3:k,2 " (512C23 + e'^' C12S23S13) (c/iH,(i) /3-k,i + e'' ^i^2(i) <i) + 

+ 523C13 {U^iit) PL^^, ~ V^,{t) al^lg) , (16) 

/3:k,r(0 = C23C13 /3!:k,3 " (^12^23 + e"^^ S12C23S13) (u^^{t) PL^^^ + ^''^23(0 

+ (si2S23 - e-^' C12C23S13) {U^t) pL^^, + e'^l/i'^3(t) a;,^i) . (17) 

These operators satisfy canonical (anti)commutation relations at equal times. The main difference with respect to 
their "naive" quantum- mechanical counterparts is in the anomalous terms proportional to the Vij factors. In fact, 
Ufj and are Bogoliubov coefficients defined as: 



t^fc.i + toA 2 / Wfe J + TOj V 



V 2uj 



k.j 



{cjk.i + mi){ujk,j + TJij) 

Ikl Ikl 



(tJ/cj+TOj) (w/c^i+TOj) 



cos(eH) 



sin(C»'^-: 



(19) 
(20) 



+ = 1 (21) 

where i, j = 1, 2, 3 and j > i. The following identities hold: 

V^Mmfit) + U^,*m^,{t) = U^,(t) , V,%{t)U^,*{t) - U^,*{t)V,%{t) = -VtAt) (22) 

u^2m^3it)-v^^2it)v,%{t) = uut) , u^3mh{t) + u^2*m2Mt) = vMt) (23) 

v.fmMt) + u^2*muit) - ^^2*^3(0 , vt2m^,it) ~ u^2mt3it) = -v,\{t) , (24) 

^5-3 = ^12 + & , e!5- = arctan {\V^\ / |C/^.|) . (25) 



As already observed in Ref. [y we remark that, in contrast with the case of two flavor mixing, the condensation 
densities are now different for particles of different masses: 
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j^i = f{m\N^mt))f = fim\N^;'\oit))f = sIAs 1^21' + inli' , (26) 

^2 = fm)\N^mt))f- fm\N^f\0it))f = \-s,,c,^ + e^'c,2S23Si3\' \V^2\' + shcl,\V^,\', (27) 
A/? = /(0(t)|A^„YlOW)/= /(OWI<riOW)/ = 1-^12^23 + e*'si2C23Si3f|F^ (28) 



nil 


m2 


ma 


6*12 


dl3 


6'23 


(5 


1 


200 


3000 


7r/4 


7r/4 


7r/4 


7r/4 



Table 1: The values of masses and mixing angles used for plots 
We plot in Fig.(l) the condensation densities for sample values of parameters as given in the table abov( 




1 10 100 1000 10000 

k 



Figure 1: Plot of the condensation densities in function of |k| for the values of parameters as in Tab.(l). 



III. THE PARAMETERIZATIONS OF THE THREE FLAVOR MIXING MATRIX 



In §11 we have studied the generator of the mixing matrix U of Eq.(g). However, this matrix is only one of the 
various forms in which a 3 x 3 unitary matrix can be parameterized. Indeed, the generator Eq.(^ can be used for 
generating such alternative parameterizations. To see this, let us first define in a more general way the generators 
Gij including phases for all of them: 



Gi2(t) = exp 6*12^12 ((5i2,t) 

G23(0 = exp 6*23 ^23 ((523, i) 

G'i3(t) = exp 6'i3Li3((5i3,i) 



-^23 ((523 I i) 



d-^x 



4{x)iy3{x)e-'^^'> ~ vl{x)iy2{x)e'^^^ 



Li3{Si3,t) = / d^x Ll{x)i^3{x)e-'''' - 4{x)i^i{x)e'''^ 



(29) 
(30) 
(31) 



■^Here and in the following plots, we use the same (energy) units for the values of masses and momentum. 
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Six different matrices can be obtained by permuting the order of the Gij (useful relations are listed in Appendix 
A) in Eq.(|). We obtain: 

12 

C12C13 si2Ci3e-**i^ sise"^'^ 

Sl2C23e**^^ -S23Sl3Cl2e'(*i^-^^^) Ci2C23-S23Sl3Sl2e-'(*^^-^i^'+*i=) S23Cl3e-*^=^ | (32) 
Sl3Cl2C23e**^' + Sl2S23e'^^^^+^^^^ -C23Si3Si2e'^^^^-^^-''> - S23Ci2e'^^^ C23C13 



5126-**!^ Si3Ci2e-**i3 

- S23Si3e*(*"-''=^^) C12C23 -si3C23Si2e'(*i= + S23Cl3e-^*^^ | (33) 

Sl2S23Cl3e'^^^^+^='^ -Cl2S23e'''^^ C23C13 + Si2Sl3S23e*('^^^+''=^-^") 

G3 = G13G23G12 

'C12C13 + Sl3S23Sl2e*(''i^-*i^+'^^-^) SuCiaC^'^'^ - Si3S23Cl2e'('^^^-''"^ Sl3C23e-''*^^ 

U3 = I -Sl2C23e'*l^ C12C23 5236"**^^ | (34) 

^Ci3S23Sl2e'(''^^+'^i^) - Si3Ci2e'^^' -Cl3S23Cl2e'^'' - Si2Si3e*(*'"-*^=) C23C13 



Sl2Cl3C23e-**i^ - Si3S23e^(''=^-*") Sl2S23Cl3e-^(*i^+*^^) + 5130236-^^1^ 

iS2 





C12C23 S23Ci2e '"^"^ (35) 




Cl2Sl3e'"^^ -0135236*'^^^ - Si2C23Sl3e'^^''''^'^'> C23C13 " Sl2S23Sl3e-'('^i^+*^^-'^") 



Sl2C23e-**l^ - Si3Ci2S23e-*^*l^-'^^^) 5130120236"*'^" + 5125236-^(^1^+*^^)' 

612623+5125235136^(^1^-*"+'^'^='' 523ei26-*'^^3 - 5i2623Sl36*(*"-''") | (36) 
-6135236"^^=^ 623613 

Gg = G12G23G13 

/ 612613 - Si2523Sl3e-*('*"+'*^3-*") 5i26236-*''i= 6125136-*'*" + 5i2S236l36-''('^i^+'^^^) \ 

Ue = -6125235136*^'^"-*^=') - 5126136**^^ 612623 6125236136-**^=' - 5i2Si36*(*"-*") (37) 

V -623S136**" -5236**^=5 623613 / 

The above matrices are generated for a particular set of initial conditions, namely for those of Eq.(|^). The freedom 
in the choice of the initial conditions reflects into the possibility of obtaining other unitary matrices from the above 
ones by permuting rows and columns and by multiplying row or columns for a phase factor. 



We thus can easily recover all the existing parameterizations of the CKM matrix |£1 24 - 29[| : 

- the Maiani parameterization (2jj2^ is obtained from I4i by setting 6*12 9, 613 /3, 6*23 7, S12 — > 0, S13 0, 
i523 -<5; J— irri 

- the Chau-Keung parameterization I2jj2^ is recovered from Ui by setting S12 and S23 0; 

- the Kobayshi-Maskawa pl| , p^ is recovered from by setting 6*12 6*2, 6'i3 9i, 623 O3, S12 —S, 613 
and S23 0, 9i ^ ^tt — 6i, with i = 1, 2, 3, and multiply the last column for (—1); 

- the Anselm parameterization [^,^ is obtained bomUi by setting 6*12 6*13, then 612 0, S13 — > 0, ^12 7r + 6'i2, 
^13 ^ TT — 6*13, 6*23 ^ |7r + 6*23, exchanging second and third column and multiplying the last row for (—1). 

From the above analysis it is clear that a number of new parameterizations of the mixing matrix can be generated 
and that a clear physical meaning can be attached to each of them, by considering the order in which the generators 
Gij act and the initial conditions used for getting that particular matrix. 



IV. CURRENTS AND CHARGES FOR THREE FLAVOR FERMION MIXING 

In this Section we study the currents associated to the Lagrangians Eq.(|^) and (|l|). To this end, let us consider the 
transformations acting on the triplet of free fields with different masses in the line of Rcf. B. 
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£ is invariant under global U{1) phase transformations of the type = e*" as a result, we have the 

conservation of the Noether charge Q = J (fix I^{x) (with /^(x) = ^m{x) 7^ ^'^(a;)) which is indeed the total charge 
of the system (i.e. the total lepton number). 

Consider then the SU{3) global transformations acting on "iim- 

^'^{x) = e'"^-^^ *,„(a;) , j = 1, 2, 8. (38) 

with Oj real constants, Fj = being the generators of SU{3) and Xj the Gell-Mann matrices pOt . 
The Lagrangian is not generally invariant under ( p8[ ) and we obtain, by use of the equations of motion, 

SC{x) = iaj^^{x)[Fj,Md]'ifm{x) = -aj d^,J^ .j{x) 

Jn^A^) = ^Mj^Fj'^^ix) , J = 1,2,.. .,8. (39) 
It is useful to list explicitly the eight currents (we suppress spacetime dependence): 



•^™,5 = [i^i7''i^3 - y^l'^yi] , J,';,,6 = ^ [i^27^'^3 + i^37^i'2] (42) 



^™,7 = -^[i^27^i^3 - ^^37''i'2] , J^,8 = ^[^^17^^^! + i^27^'^2 - 2P37''^3] ■ (43) 

The related charges Qm.jit) = / d'^x J^ ,j{x), satisfy the su{3) algebra [Qm,j{t),Qm.k{t)] = i fjki Qm,i{t)- Note that 
only two of the above charges are time- independent, namely Qm.,3 and Qm.,8- We can thus define the combinations: 

Qi = + Q,„,3 + -^Qm,8, (44) 

Q2 = ^Q ~ Q,n,3 + ;^*5m,8, (45) 
Q3 = ^Q - ^Qm,8, (46) 



, z = l,2,3. (47) 



These are nothing but the Noether charges associated with the non-interacting fields vi , V2 and : in the absence of 
mixing, they are the flavor charges, separately conserved for each generation. 

As already observed in §11, in the case when CP is conserved (5 = 0), the mixing generator Eq.(||) is an element of 
the SU (3) group and can be expressed in terms of the above charges as: 

Ge{t)\^_^ = gi2e23Q,„,7(t) g»2ei3Q,„,5(t) g»2ei2Q™,2(t) (-43-) 

Following Ref. Q|, we can now perform the SU{2i) transformations on the flavor triplet 4"/ and obtain another set 
of currents for the flavor fields: 

*'/x) = e'"^^^*/(a;) , j = l,2,...,8, (49) 

which leads to 

5C{x) = za,l'/(x)[F„M]*/(x) = -a^d^J'^^x), 

J';Ax) = ^'f{x)^^F,^>f{x) , j = l,2,...,8. (50) 
Alternatively, the same currents can be obtained by applying on the ^ {x) the mixing generator Eq. (|^) : 
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JfA^) = Ge\t)J':^^{x)Gg{t) , j = l,2,...,8. (51) 

The related charges Qf,j{t) = / S'xJ'j ^{x) still close the su(3) algebra. Due to the off-diagonal (mixing) terms in 
the mass matrix M, Qf_z{t) and Qf^s{t) are time-dependent. This implies an exchange of charge between Vg, and 
Pr, resulting in the flavor oscillations. 

In accordance with Eq. ( ^4[ ) , we define the flavor charges for mixed fields as 

Qr{t) ^Iq- -^Qf^^i*)- (54) 
with Qe{t) + Quit) + Qrit) — Q- These charges have a simple expression in terms of the flavor ladder operators: 
Q.it) -Y.h'^ (<.(i)«L,.(i) - /3:U(0/3:k,.(i)) , ^ = e, T , (55) 

because of the connection with the Noether charges of Eq.(^7|) via the mixing generator: Qa{t) — Gg^{t)QiGe{t). 

Notice also that the operator AQa-{t) = Qa{t) — Qi with (a, i) = (e, 1), (/i, 2), (r, 3) , describes how much the mixing 
violates the (lepton) charge conservation for a given generation. 

Let us now come back to the algebra of the currents including CP violating phases. To this end, we consider a 
generalization of the Gell-Mann matrices (we use a tilde for denoting the modified quantities including phases): 

e'^^ 0\ / -ie'*^ 0\ /l 

e-''*2 , A2 = ie-"^'- , A3 = -1 | , A4 
0/ \ 0/ \0 

-ie^^^A /O \ /O 

A5=l0 ,A6= e'^' , A7 = -ie"^'^ I Ag = ^ I 1 I . (56) 




vie 



0/ \0 e-'*^ / \0 ie-'^-' 
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These are normalized as the Gell-Mann matrices: tr{XjXk) = 26jk ■ We define as usual the algebraic generators: 



^. = ^A, , j = l,..,8 (57) 



The above generators Eq.(p7|) do not close the su(3) algebra, unless the condition 62+65 + 67 = is imposed 
(cf. Eq.(^8|) below), i.e. if one of the three phases is fixed in terms of the remaining two. Such a request is clearly 
incompatible with the parameterizations of the mixing matrices of §11 and §111 (cf. e.g. the discussion after Eq. (^7|) ; 
we have the correspondence {62, 65, 67} {612, 613, 623})- For the matrix (@), as already observed, it implies 6 — 0. 

The Fj satisfy a deformed su(3) algebra with deformed commutation relations given by: 



[F2,F5] 


^\F7e 


-iA{Fa-V3F8) . 


[^2,^7] = - 


-he 
2 ^ 


-iA{F3 + V3Fs) . 


[F5,F7] 




[Fi,F4] 


i ~ 
^-^F.e- 


-i A(F3-\/3Fs) . 


[Fl.Fj] 


i ~ 


'lAiFs + VSFa) . 


[Fi,F7] 


"2^ 


[Fi,F5] 


2 ^ 




; [Fi,Fe] = 




-iAiFa + VSFg) . 


[F5,Fe] 


- 2^ie 


[F2,Fi] 




-iA{Fa~V3F8) . 


[F2,F6] = - 


\he 


-iA{F3 + V3F8) . 


[F4:,Fg] 


2^ 



(58) 

where A = 62 + 65 + 6j. The other commutators are the usual su{3) ones. For A = 0, the su{3) algebra is recovered. 
It is useful to look at the deformed algebra in terms of the raising and lowering operators, defined as [ pO| : 

f± = Fi±iF2 , U± = Fe±iF7 , V± = F4 ± (59) 
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We also define: 

Then the only deformed commutators are the following ones: 



F3 = i {V3Fs 



2iAU3 



[f+,U+] = v+ 



[U+,V.] = f_e 



2iAT, 



(60) 



(61) 



V. NEUTRINO OSCILLATIONS 



The oscillation formulas are obtained by taking expectation values of the above charges on the (flavor) neutrino 



state. Consider for example an initial electron neutrino state defined as \ve) = Q!k^e('-')l^)/ (^°^ * discussion on the 



correct definition of fiavor states see Refs. PfllTl]]). Working in the Heisenberg picture, we obtain 

Ql^t) EE (^pig.(i)i^p) - /(o|g.(t)|o)/ = |{"L.W,<p(o)}f + |{/3:l.(i),<p(o)}f , (62) 
QiAt) ^ {^p\Q^{t)K) - /(0|Q.(t)|0)/ = - |{/3L.(0,/3k!p(0)}f - |{"-l.W'/3k!p(0)}f , (63) 

where |0)/ = |0(0))/. Overall charge conservation is obviously ensured at any time: Qk,e(i) + Qk,/j(0 + 2k,T(^) = 1- 
We remark that the expectation value of Qa cannot be taken on vectors of the Fock space built on |0)m, as shown in 
Refs. Also we observe that /(0|Qcr(t)|0)/ 7^ 0, in contrast with the two flavor case P;p[- We introduce the 

following notation: 



A' 



i^k.: 



Then the oscillation (in time) formulas for the flavor charges, on an initial electron neutrino state, follow as: 
Qk,e(^) = 1 - sin2(20i2)cos*0i3 \\U^2? sin^ (A^ai) + \V^2? sin^ [n^^t] 



sin^ (26*13) cos^ O12 



C/l'aP sin^ (A^^^i) + |T^i3p sin^ (^^^'3* 



^13 



13'' 



Sin2(26'i3)sin2 6'i2 \U^3\'^ sin^ (A^3i) + IF23P sin^ {9^^t) 



Qw(i) = 2J,, sin(2A^20 - \V^\' sin(2f]^2t) + - |l^i'^3r) sin(2A^3i) 



(64) 



+ (1^21' - |Fi''3r)sin(2f7^3t) - |[/f3|^ sin(2A^3i) + \V^^\' s\n{2a%t) 



cos ^13 sm (^13 



cos Wi3SmC'i3 



f/^3psin2(A^3i)+|^i'^3p sin2(l7^3t) 



COS 6 sin(20i2) sin(26'23) + 4 cos^ 6x2 sin ^13 sin^ 6*23 

cos(5sin(26li2) sin(26'23) - 4sin2 612 sin 6*13 sin^ 6123] \U^:i\^ sir? {iS^^^t) + \V^^\^ sir? {^^t 
cos^ 6^13 sin(26'i2) (cos^ 6*23 - sin^ 6*23 sin^ 6*13) sin(26'i2) + cos(5cos(26li2) sin 6*13 sin(26'23) x 
U^2 1 ' sin^ ( A^2i) + I Vt2 1 ' sin^ {^\^t) ] , (65) 

Qf{t) ^ -2J,, \\U^,\' sin(2A^2i) - IV.W' sm{2n'l,t) + {\U^,\' - \V,%n sm{2A^,t) 



k |2 



Fi^3|2)sin(2f]^3i) - \U^,\-' sin(2Af30 + \V,%{' sm{2n'^,t) 



rk |2 



k |2 



cos t/i3 sm 



COS c/13 sm6'i3 



\U^,\' sin' (A^3t) 



U^,\' sin^(A?3i) + |K,'^|^ sin^(r!-i 



^231 



^23 



cos(5sin(26'i2) sin(2023) — 4cos^ 612 sin0i3 cos^ 6'23 
cos(5sin(26'i2) sin(26'23) + 4sin^ 6*12 sin 6*13 cos^ 6*23 
+ cos^ 6*13 sin(20i2) (sin^ ^23 — sin^ ^13 cos^ 6*23) sin(2f?i2) — cos (5 cos(20i2) sin6'i3 sin(2023) 
\U^2\' Sin' (A^2t) + |Vi^2l' sin^ {n^,t) 



v,\\' sin^ {n^,,t) 



k |2 



231 



'23'' 



(66) 
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where we used the relations Eqs.(g2|)-(g5|). We also introduced the Jarlskog factor J^p defined as I^Oj 

= Im{u,aUjf3U*pU*^), (67) 
where the Uij are the elements of mixing matrix lA and i ^ a ^ j3. In the parameterization Eq.(|^), J^p is given by 



Jcp = \ sin5 sin(26'i2) sin(26'i3) cos 6*13 sin(26'23)- 



(68) 



Evidently, J^^ vanishes if 6ij = 0, 7r/2 and/or 5 — Q,tt: all CP-violating effects are proportional to it. 

The above oscillation formulas are exact. The differences with respect to the usual formulas for neutrino oscillations 
are in the energy dependence of the amplitudes and in the additional oscillating terms. For |k| ^ y^rriim^, we have 
\U^\'^ — > 1 and \V^\'^ and the traditional (Pontecorvo) oscillation formulas are approximately recovered. Indeed, 
for sufficiently small time arguments, a correction to the Pontecorvo formula is present even in the relativistic limit. 

In Appendix B the oscillation formulas for the fiavor charges on an initial electron anti-neutrino state are given. 

We plot in Figs. (2) and(4) the QFT oscillation formulas Qke(^) ^^'^ 2k/j(^) ™ function of time, and in Figs. (3) 
and (5) the corresponding Pontecorvo oscillation formulas P^_,f^{t) and P^^(t). The time scale is in T12 units, where 
2^12 = '^/^^2 is, for the values of parameters of Tab.(l), the largest oscillation period. 



. 5 




Figure 2: Plot of QFT oscillation formula: ^{t) in function of time for A; = 55 and parameters as in Tab.(l). 
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Figure 5: Plot of QM oscillation formula: P^_^^{t) in function of time for fc = 55 and parameters as in Tab.(l). 



VI. CP AND T VIOLATIONS IN NEUTRINO OSCILLATIONS 



In this Section we consider the oscillation induced CP and T violation in the context of the present QFT framework. 
Let us first briefly recall the situation in QM^: there, the CP asymmetry between the probabilities of two conjugate 
neutrino transitions, due to CPT invariance and unitarity of the mixing matrix, is given as pq] 



where a, p = e, fi, t. The T violating asymmetry can be obtained in similar way as [p5| 



(69) 



(70) 



The relationship A^J(i) = A^'^(t) is a consequence of CPT invariance. 

The corresponding quantities in QFT have to be defined in the framework of the previous Section, i.e. as expectation 
values of the flavor charges on states belonging to the flavor Hilbert space. We thus have for the CP violation: 



(71) 



{"L,.(0,<p(0)}f + |{/3:U(0,<p(0)}f - \{aXMf^llm}f |{/3k,.(0,<p(0)}f • (72) 



We have 



p,cr = e,^,T, 



(73) 



which follows from the fact that J^aQr^i^) ~ Q ^^'^ i'^plQWp) = 1 ^.nd (t'plQjt'p) = —1. 
We can calculate the CP asymmetry Eq.([7l|) for a specific case, namely for the transition 



k |2 



k |2 

2l 



k |2^ 
13 1 / 



sin(2A^3t) 



-{\Vh\' - \V,\nsm{2n^,t) - \U^,\' sin(2A^3t) + \V,\\' sin(2f7^3t) 



z/p. We obtain 



(74) 



We use here the "hat" for QM quantities. For notational simplicity, we also suppress momentum indices where unnecessary. 
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and A^'^(<) = —A'^^{t). As already observed for oscillatfon formulas, high-frequency oscillating terms and Bogoliubov 
coefficients in the oscillation amplitudes appear in Eq. ( [74|) as a QFT correction to the QM formula. 

The definition of the QFT analogue of the T- violating quantity Eq.(70) is more delicate. Indeed, defining as 
A^'' = Q^(i) — Qe (0 does not seem to work, since we obtain A^'' — A^^ ^ in contrast with CPT conservation. 

A more consistent definition of the time-reversal violation in QFT is then: 



k.cr 



p,a ^ e, ^, 



(75) 



With such definition, the equality A^'"(i) = MZit) follows from Q£,<,(-t) = -Q^,<,(t)- 

We plot in Fig. (6) the CP asymmetry Eq.(|7J) for sample values of the parameters as in Tab.(l). 
corresponding standard QM quantity is plotted for the same values of parameters. 



In Fig. (7) the 



0.5 - 




-0.5 



Figure 6: Plot of the QFT CP asymmetry A'^{t), in function of time for fc = 55 and parameters as in Tab.(l). 



0.5 - 



-0.5 




Figure 7: Plot of the QM CP asymmetry A^'^(t), in function of time for fc = 55 and parameters as in Tab.(l). 
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VII. CONCLUSIONS 



In this paper, we have discussed the mixing of (Dirac) fermionic fields in Quantum Field Theory for the case of 
three flavors with CP violation. Wc constructed the flavor Hilbcrt space and studied the currents and charges for 
mixed fields (neutrinos). The algebraic structure associated with the mixing for the case of three generation turned 
out to be that of a deformed su{3) algebra, when a CP violating phase is present. 

We have then derived all the known parameterization of the three-flavor mixing matrix and a number of new ones. 
We have shown that these parameterizations actually reflect the group theoretical structure of the generator of the 
mixing transformations. 

By use of the flavor Hilbert space, we have calculated the exact QFT oscillation formulas, a generalization of the 
usual QM Pontecorvo formulas. The comparison between the exact oscillation formulas and the usual ones has been 
explicitly exhibited for sample values of the neutrino masses and mixings. CP and T violation induced by neutrino 
oscillations have also been discussed. 

We remark that the corrections introduced by the present formalism to the usual Pontecorvo formulas are in 
principle experimentally testable. The fact that these corrections may be quantitatively below the experimental 
accuracy reachable at the present state of the art in the detection of the neutrino oscillations, does not justify 
neglecting them in the analysis of the particle mixing and oscillation mechanism. The exact oscillation formulas here 
derived are the result of a mathematically consistent analysis which cannot be ignored in a correct treatment of the 
field mixing phenomenon. As we have seen above, our formalism accounts for all the known parameterizations of the 
mixing matrix and explains their origin and their reciprocal relations, thus unifying the phenomenological proposals 
scattered in the literature where such parameterizations have been presented. Moreover, our formalism clearly points 
to the truly non-perturbative character of the particle mixing phenomenon. A lot of Physics must be there waiting 
to be discovered. 



ACKNOWLEDGEMENTS 



The present research has been carried out in the framework of the ESF network COSLAB and we acknowledge 
partial financial support also by INFN, INFM and EPSRC. 



APPENDIX A: ANTI-NEUTRINO OSCILLATION FORMULAS 



If we consider an initial electron anti- neutrino state defined as Ive) = /^k^e(^)l^)/' ^® obtain the anti- neutrino 
oscillation formulas as 



(76) 



Qijt) = 2J,, \U^^\' sin(2A^2i) - K2V ^H'^^l^t) + {\U^2V ' \V^\')sm{2i^t^t) 



|yi''3nsin(2f7^3t) - |;7^3p sin(2A^3t) + \V^^\' sm{2n^,t) 



cos (^i3sin6'i3 



U^,\'sin' (A^3t) 



+ cos ^13 sin 



cos 6 sin(20i2) sin(2^?23) + 4 cos^ 612 sin Oi^ sin^ 623 

cos5sin(26»i2)sin(26'23)-4sin2 6li2sin6'i3sin2 6l23] IC/^sP sin^ (A^3t) + 1^/23 P sin^ (f^^3t) 
— cos^ ^13 sin(2^i2) (cos^ 623 — sin^ 623 sin^ ^13) sin(2^i2) + cos^ cos(2^i2) sin ^13 sin(2^23 
\U^2\^ sin^ (A^^t) + \V^2\^ sin2 [n^^t) 



(77) 
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cos t/i3 sm 



COS t/i3Sin6'i3 



cos(5sin(26'i2) sin(2023) — 4cos^ 6*12 sin^ia cos^ 6'23 



iC/i 



k |2 
13 1 



sin^ (A^,t) + \V, 



k |2 

13 1 



sm 



cos(5sin(26'i2) sin(2023) + 4sin^ 612 sin^is cos^ 



23 



cos' ^13 sin(20i2) (sin^ ^23 - sin^ 613 cos^ ^23) sin(20i2) - cos 5 cos(20i2) sin ^13 sin(2023) 



\U^2\' sin' (A^2t) 



|Vi^2l' sin' (n^^t) 



(78) 



APPENDIX B: USEFUL FORMULAS FOR THE GENERATION OF THE MIXING MATRIX 



In deriving the Ui mixing matrices of §11 and §111, we use the following relationships 

K"(a:),Li2] = u^{x)e-''-^ , K"(a:),i23] = , K(x),Li3] = ,^i{x)e-''^^ , 

K(x),Li2] = -i^n^y''' , K(^),^23] = i^iix)e-^'^^ , K(x),Li3] = , 



and 



G2-3Hi);^r(a;)G23(i) 

Gr3^W;.n2^)Gi3W 



z.i"(a;)ci3+z/3"(a;)e-'^i^si3 , 
i'r(a;)ci2 + i^2"(a^)e"''^''si2 , 



(79) 
(80) 
(81) 



(82) 
(83) 
(84) 



G2-/(t)z.2"(x)G23(t) 

Gr/(t)z.2"(a:)Gi3(t) 



V2{x)c2s + V3{x)e 
v^{x) , 

V2{x)ci2 - v'^{x)e 



-i&2 



S23 



S12 , 



(85) 
(86) 
(87) 



G2-3i(t)z.3"(a;)G23(i) 
Gr3'(t);.3"(a;)Gi3(i) 
Gr2'(t);.3"(a;)Gi2(i) 



1^3 (a;)c23 - V2{x)e'^'-S23 , 
y^{x)c^3-vi{x)e'''s,3. 
1^3 



(89) 
(90) 



APPENDIX C: ARBITRARY MASS PARAMETERIZATION AND PHYSICAL QUANTITIES 

In RefjBQ it was noticed that expanding the flavor fields in the same basis as the (free) fields with definite masses 
(cf. Eq.(|ll[)) is actually a special choice, and that a more general possibility exists. In other words, in the expansion 
Eg. (pi]) one could use eigenfunctions with arbitrary masses ^a, and therefore not necessarily the same as the masses 
which appear in the Lagrangian. On this basis, the authors of Ref. have generalized the Blasone-Vitiello (BV) 
formalism by writing the flavor fields as 



(91) 



where u„ and are the helicity eigenfunctions with mass /io- . We denote by a tilde the generalized flavor operators 
introduced in Ref. ||^,^ in order to distinguish them from the ones in BV formalism Eg. (pi]). The expansion Eg.(|9l|) 
is more general than the one in Eg.([l^) since the latter corresponds to the particular choice /ie = toi, /i^ = ?Ti2, 
/j,T- = TO3. Of course, the flavor fields in Eg. (|9l|) and Eg.(|l^) are the same fields. The relation, given in Ref. 
between the general flavor operators and the BV ones is 
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with (cr,j) = (e, 1), (^,2), (t,3), £^ ^ = (x^ - Xj)/2 and cotXa = |k|/M^, coix) = |k|/TOj. For = m^, one has 

As aheady noticed in Ref. g, the flavor charge operators are the Casimir operators for the Bogohubov transfor- 
mation (|9^), i.e. they are free from arbitrary mass parameters : Qa{t) = Qait)- This is obvious also from the fact 
that they can be expressed in terms of flavor fields (see Ref. |Q ) . 

Physical quantities should not carry any dependence on the ^a-'- in the two-flavor case, it has been shown |^ that the 
expectation values of the flavor charges on the neutrino states are free from the arbitrariness. For three generations, 
the question is more subtle due to the presence of the CP violating phase. Indeed, in Ref. it has been found that 
the corresponding generalized quantities depend on the arbitrary mass parameters. 

In order to understand better the nature of such a dependence, we consider the identity: 

im^m) = (^1 j(o) Q^{t) j-\om = {mAtm + (^i [^w, Qam j-\om) . (94) 

valid on any vector of the flavor Hilbert space (at t — 0). From the explicit expression for J(0) we see that the 
commutator [J(0), Qa-it)] vanishes for /ip = irij , {p,j) = (e, 1), (/i, 2), (r, 3). 
It is thus tempting to define the (effective) physical flavor charges as: 

Qf'y^t) = Q,{t) - J-i(O) [J(0), Q,(t)] - J-i(O) Q,{t) J(0), (95) 

such that for example: 

{^p\Ql''''{t)\i^p) = {^p\Q.{t)K)- (96) 

It is clear that the operator Q^y^{t) does depend on the arbitrary mass parameters and this dependence is such 
to compensate the one arising from the flavor states. The choice of physical quantities (flavor observables) as those 
not depending on the arbitrary mass parameters is here adopted, although different possibilities are explored by other 
authors, see Refs. [^jTUTl. 
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